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In a previous work we showed that massive test particles exhibit a non- 
geodesic acceleration in a modified theory of gravity obtained by a non- 
commutative deformation of General Relativity (so-called Matrix Grav- 
ity). We propose that this non-geodesic acceleration might be the origin of 
the anomalous acceleration experienced by the Pioneer 10 and Pioneer 1 1 
spacecrafts. 
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1 Introduction 

The Pioneer anomaly has been studied by many authors (see [HI |2l [TTl [121 El 
[T3l and the references in these papers) and it has a pretty strong experimental 
status [[91- It exhibits itself in an anomalous acceleration of the Pioneer 10 and 1 1 
spacecrafts in the range of distances between 20AU and 50AU (~ lO^'^cm) from 
the Sun. The acceleration is directed toward the Sun and has a magnitude of [[B 1211 

^anom « (8-74 ± 1.33) X lO'^cm/s^ . (1.1) 

In the last years there have been many attempts to explain the Pioneer anomaly by 
modifying General Relativity (see, for example, [fTTll and the references therein). 
However, there is also some evidence [I13II that it could not be explained within 
standard General Relativity since it exhibits a non-geodesic motion. That is, it 
cannot be explained by just perturbing the Schwarzschild metric of the Solar sys- 
tem. It seems, from the analysis of the trajectories, that the spacecrafts do not 
move along the geodesies of any metric. Another puzzling fact is that there is 
no measurable anomaly in the motion of the planets themselves, which violates 
the equivalence principle. In other words, the heavy objects like the planets, with 
masses greater than ~ lO^^g, do not feel any anomaly while the smaller objects, 
like the Pioneer spacecrafts, with masses of order ~ lO^g, do experience it. 

There are also some interesting numerical coincidences regarding the Pioneer 
anomaly (noticed in [10] as well). Recall that the cosmological distance, which 
can be defined either by the Hubble constant H or by the cosmological constant 
A, is of order 

ro-^--^- lO^^cm (1.2) 
H Va 

and the Compton wavelength of the proton is of order 

^ J_ ^ lO-i^cm. (1.3) 

MpC 

Now, we easily see, first of all, that there is the following numerical relation 

ri \ /ranomf^ ^^^^ 



where ranom ~ lO'^^cm is the distance at which the anomaly is observed. This 
means that 

h I He 
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Secondly, the characteristic distance determined by the value of the anomalous 
acceleration, Aanom ~ lO^^cm/sec^, is of the same order as the cosmological dis- 
tance 

r2 lO^^cm, (1.6) 

■'4 anom 

which simply means that 

A,,,^-^ Hc-^c^yfX. (1.7) 

It is very intriguing to speculate that the Pioneer effect is the result of some kind 
of interplay between the microscopic and cosmological effects at the macroscopic 
scales. 

In this paper we apply the investigation of motion of test particles in an ex- 
tended theory of gravity, called Matrix Gravity, initiated in f6l to study the anoma- 
lous acceleration of Pioneer 10 and Pioneer 11 spacecrafts. Matrix Gravity was 
proposed in a series of recent papers [I3l|4l[5l. This is a modification of the stan- 
dard General Relativity in which the metric tensor g'^'^ is replaced by a Hermitian 
NxN matrix-valued symmetric two-tensor af^^ = g^^l-irh^^ , where I is the identity 
matrix, is a matrix- valued traceless symmetric tensor, i.e. \xW^ = 0. In this 
theory the usual interpretation of gravity as Riemannian geometry is no longer 
appropriate. Instead, Matrix Gravity leads, quite naturally, to a generalized geom- 
etry, that we call Matrix Geometry, which is equivalent to a collection of Finsler 
geometries. Instead of a usual Riemannian geodesic flow, we get a system of 
Finsler flows, and, moreover, the mass of a test particle is replaced by a collection 
of mass parameters. In the commutative limit, only the total mass is observed. For 
more details and discussions see Il4l|5l|6]|. 

The dynamics of the tensor field a'^'' is described by a non-commutative Einstein- 
Hilbert action, which can be constructed either by an extension of all standard 
concepts of diff'erential geometry to the non-commutative setting dH HI or from 
the spectral invariant of a partial differential operator of non-Laplace type 

The main goal of the present paper is to apply our previous study flU of the mo- 
tion of test particles (in a simple model of matrix gravity) to the Pioneer anomaly. 
We would like to stress that this study is just a first attempt to analyze the phe- 
nomenological effects of Matrix Gravity. We do not claim that this simple model 
definitely solves the mystery of the anomaly. Our aim is just to propose another 
candidate for its origin. Only future tests and more detailed models can describe 
the Pioneer anomaly in full capacity. This work does not represent the final an- 
swer, but just a first attempt of studying this phenomenon within the framework 
of Matrix Gravity. 
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2 Anomalous Acceleration in Matrix Gravity 

In Matrix Gravity a massive particle is decsribed not by a single mass parameter 
m but rather by A'^ different mass parameters m,, so that m = ^jli ^i- I^i the 
commutative limit we only observe the total mass m. The interesting question of 
the physical origin of the parameters m, requires further study. For this reason, we 
do not assume that the m, are positive. Following [6J we consider two different 
cases. In the first case, that we call the nonuniform model, we assume that mass 
parameters are different, and in the second case, that we call the uniform model, 
we discuss what happens if they are equal to each other, that is, m, = mfN. 

The equations of motion of a test particle are derived and studied in |[6l. They 
have the form ^ 

^+rVx,i)i"/ = 0, (2.1) 

where y^Q.^(jc, x) are generalized Christoffel coefficients that are homogeneous 
functions of x of order zero, in other words, they depend on the direction of x, but 
not on its magnitude. These functions depend in a complicated way on the matrix- 
valued metric a'^^, on the velocity, xf, and, in general, on the ratios //, = m,7m. 

In the perturbation theory, when one writes a^^ = gf^^l + h'^", the generalized 
Christoffel coefficients are 

rafiix, X) = r^^(x) + (T^pix, X) , (2.2) 

where Y^ap are usual Christoffel coefficients of the metric of the metric g^y and 
ff^a/ii^, ^) is some tensor of first order in the perturbation. Now, the matrix-valued 
metric a^^ = g'''^! + h^^ satisfies the non-commutative Einstein equations derived 
in [|4l[51|7]|- As a result of non-commutative corrections even the equations for the 
metric g'^'*' are modified. This means that both the metric g^"" and the Christoffel 
symbols T^a/i are modified. More precisely, we let 

g'''' = g^;+r, (2.3) 

where is the non-perturbed Riemannian metric given by the solution of the 
standard Einstein equations without any non-commutative corrections and jS'^'" is 
the non-commutative correction. Then 

r^/j = r^.^ + a^/^ , (2.4) 

where Fqq,^ is the Christoffel symbols of the metric g'^^ and a^^is is the perturbation. 
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Thus, the equations of motion take the form 




(2.5) 



where 



non-geod 



{x, x) 



(2.6) 



is the anomalous acceleration and and 



A'l^^^{x,x) = -a^,p{x)x"^, 

on-geod(^' = -ff'apix, x)x" Xp . 



(2.7) 



(2.8) 



are the geodesic and non-geodesic parts of the anomalous acceleration. It is this 
anomalous acceleration that we are going to study in this paper. We suggest that 
this might explain the anomalous behavior of Pioneer 10 and 1 1 spacecrafts. 

The anomalous geodesic acceleration can be easily computed by expanding 
the Christoffel coefficients in the perturbation. In the first order in yS we obtain 



where the covariant derivatives and all tensor operations are performed with the 
non-perturbed metric g^^. 

The anomalous non-geodesic acceleration was derived within perturbation 
theory in the deformation parameter in We study the two cases mentioned 
above. 

Nonuniform Model. First, we study the generic case when the parameters 
are different. We define a function 



where is a covector and Ai{x, ^) are the eigenvalues of the matrix h'^^(x)^fj^y. 
Note that since trh^^ = the matrix h^^^^^y is traceless, which implies that the 
sum of its eigenvalues is equal to zero. Thus, in the uniform case, when all mass 
parameters /i, are the same, the function P(x,^) vanishes. In this case the effects 
of non-commutativity are of the second order. 

The non-geodesic acceleration was computed in H and has the form 




(2.9) 



N 




(2.10) 




(2.11) 
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where ^ 

qnx,^) = \-^P{x,^) (2.12) 

and the covariant derivatives are defined with the Riemannian metric. Thus, the 
total anomalous acceleration is 

A'^anom = \ [^"^ aifi" p + p) ' "^"(fia/S + Qafi)] . (2.13) 

Uniform Model. Now, we will simply assume that all mass parameters are 
equal, that is, nii = jj . The non-geodesic acceleration, computed in [6], has the 
form ^ 

A^non-geod = ~g,?'^ {^'^^ aS fjypa ~ vS a/jpaj X X xP , (2.14) 

where 

5^""^ = ^tr(m''^). (2.15) 
Thus, the total anomalous acceleration is 

(2.16) 

In the spherically symmetric background, in the non-relativistic limit, the ra- 
dial anomalous acceleration is given by: in the uniform model, 

A'-a„o. = |:(-^rw + ^5«°°V)), (2.17) 

and, in the nonuniform model, 

A^no. = |:(-^r(r)-i^°V)). (2.18) 

Of course, this can also be interpreted as a modification of Newton's Law fSJ. 

Here, of course, the tensor components /3°°, S and q^^ should be obtained 
as solution of the non-commutative Einstein field equations (in the perturbation 
theory). These equations are somewhat complicated. That is why, in this paper 
we consider a toy model just to get a glimpse into the phenomenon. 

We consider a simple model of 2 x 2 real symmetric commutative matrices. 
The static spherically symmetric solution of the matrix Einstein equations for this 
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model was obtained in [[6]|. In the spherical coordinates jc" = t, = r, = 6, 
x^ = (fit has the form 



with 



a^" = A, a'' = B, (2.19) 

a^^ = \l, a^^ = I , 

r2 r2 sin^ 6 



I . r„-eL\ L 



B{r) = 1--Ar^-^ I + -T, 



A(r) = ^(r)I + (A(r)T , (2.20) 
where 1 = 1^^ nl'^^^ cosmological constant, = 2GM, M 



1/' \ 1 ^ 

is the mass of the Sun, 6 and L are some integration parameters, (the parameter 
6 should not be confused with the angle variable 9 above), and (p(r) and ifr{r) are 
some functions (the function (p{r) should not be confused with the angle variable 
(fi above). The functions (p{r) and tff(r) can be parametrized by 

m = -Mr) + 2ef2(r), ipir) = efi{r) + (1 + e')f2{r) . (2.21) 

By introducing the following function 

u{r) = 1 - ^Ar^ - ^ , (2.22) 
3 r 

we obtain, for /i(r) and /2(r) in (|2.21|) . the expressions 

u(r) 

m = , (2.23) 



fiir) = -r . (2.24) 

[u{r) + {6 +l)^\[u{r) + {6-1)^1 

Notice that when the parameters 6 and L vanish the functions (p{r) and i/r(r) 
give nothing but the standard Schwarzschild solution with the cosmological con- 
stant CSl. We take the standard Schwarzschild solution (without the cosmological 
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constant) as the non-perturbed metric . The parameters 6, L and A will be con- 
sidered as perturbation. This means that the relevant components have the form 



1-^' 



(2.25) 
(2.26) 



and, therefore, 

SO«°o = ^V)- (2.27) 

In this simple 2x2 model the tensor q^^^ is given by q^^'^ = |tr (/z^^^r), f6l, where 
y = Hi -1^2, and therefore, 

q'' = r^(r). (2.28) 

As we already mentioned above, in the non-relativistic limit the only essential 
component of the anomalous acceleration is the radial one A''anom- By using the 
equations above we obtain: in the uniform model. 



A'' 

anom 



]_d_ 



_1_ + /j(^) _ 2ef2{r) + (^/i(r) + (1 + e^)f2(.r)^ 



, (2.29) 



and in the non-uniform model. 



anom 



ld_ 



1 



1 



j- + {ye-\)f,{r) + y{\+eYf2{r) 



(2.30) 



We would like to emphasize at this point that the perturbation theory is only 
valid for small corrections. Obviously, when the corrections become large one 
needs to consider the exact equations of motion. 



3 Pioneer Anomaly 

We have two free parameters in our model, 6 and L (and y in the non-uniform 
model). We estimate these parameters to match the value of the observed anoma- 
lous acceleration of the Pioneer spacecrafts. 

First of all, we recall the observed value of the cosmological constant A « 
2.5 • lO'^^cm"^; therefore, tq ^ |A|~'^^ = 6.3 • lO^^cm, and the gravitational ra- 
dius of the Sun Vg « 1.5 • lO^cm. The relevant scale of the Pioneer anomaly 
is /'anom ~ 10^"* ^ lO'^cm, therefore, we can restrict our analysis to the range 
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rg « r « rQ. The values of the dimensionless parameters are j ~ 10 ^, 
— ~ 10^'^, and — ~ 10"^^. We also remind that the value of the anomalous 

'•o '"0 

acceleration is A^^^^^ « 8.7 • 10~^cm/s^ . We should stress that our analysis only 
applies to the range of distances relevant for the study of the Pioneer anomaly. 
Therefore, strictly speaking, from a formal point of view, one cannot extrapolate 
our equations beyond this interval. Since the parameters ^ and ^ are negligi- 
bly small (compared to 1) they can be omitted. 

By using the eqs. (12.291 ) and (12.301 ). and by defining p = (1 + 6^)L - dr^, we 
obtain [6] (in the usual units, c being the speed of light) for « r « ro: in the 
uniform model, 

(3.1) 



and in the non-uniform model. 



(p + 29rg 2 
3' 



Knor. = ^y\—r^--(^^r\. (3.2) 



Uniform Model. First, we restrict to the case of vanishing cosmological con- 
stant. Then the function (13.11 ) takes the form 



A'"anom(r) = — T ( ^ + " ) 2 " (^.3) 



It has an extremum if the signs of 6 and p are different, which occurs at = 
and is equal to 

It ^ 

Now, we assume that ~ ranom ~ 10^"^ cm ariA A''^^^J(r^) ~ -10~*^cm/sec^ to 
estimate the parameters 



p~10'cm, e--\0~\ (3.5) 

If we leave the cosmological constant there is another range of parameters that 
should be investigated. Namely, when the term becomes comparable with 

the term ^. In this case the anomalous acceleration can be written, by dropping 
negligible terms, as 
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2 

We note that the term gives the right magnitude of the anomalous acceler- 
ation. If we assume that the two terms in the parentheses are comparable at the 
characteristic length ranom and are of order 1 , then we get an estimate 

- ''^'""■"-^ and e-i—Y , (3.7) 



'o 



and, therefore, 



p ~ 10~^cm and 6* ~ 10\ (3.8) 

Nonuniform Model. In the non-uniform model we have an additional parame- 
ter y. The function has an extremum at 



Now, we assume that 





r* = 1 


anom 




P 




e 


~ 3^ 



3P^V 



(3.9) 



^-13. 



(3.10) 



Further, by assuming A' anom('"*) 10 ^cm/sec^ and using the eq. we estimate 

the parameter y 

r-lO^^ (3.11) 

It is interesting to notice that, in this case, pl6 has the same order of magnitude 
of the Compton wavelength of the proton. Moreover, by using (13.101) . we confirm 
the coincidence (11.51) mentioned in the introduction. This is very intriguing; it 
allows one to speculate that the anomalous acceleration could be a result of an 
interplay between the microscopic and macroscopic worlds, in other words, the 
Pioneer anomaly could be a quantum effect. 



4 Conclusions 

In this paper we applied the kinematics of test particles {6\ in Matrix Gravity PUIS]! 
to the study of Pioneer anomaly. Matrix Gravity is interpreted in terms of Matrix 
Geometry, a generalized geometry which is equivalent to a collection of Finsler 
geometries, rather than Riemannian geometry. This new feature of our theory 
leads to an interesting and completely new phenomenon of splitting of Rieman- 
nian geodesies to a collection of Finsler geodesies. More precisely, instead of one 
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Riemannian metric we have different Finsler metrics and different mass parame- 
ters which describe the tendency to follow a particular trajectory determined by a 
particular Finsler metric. The interesting result is that test particles in our theory 
exhibit a non-geodesic motion which can be interpreted in terms of an anoma- 
lous acceleration. This new feature led us to apply these results for studying the 
anomalous acceleration of the Pioneer spacecrafts. 

We considered two models: a uniform one, in which a particle is described by 
a single mass parameter, and a non-uniform one, in which a particle is described 
by multiple mass parameters. The choice of one model over the other should be 
dictated by physical reasons. The interesting question of whether the matter is 
described by only one mass parameter or more than one mass parameters requires 
further study. If the Pioneer anomaly is a new physical phenomenon we have to 
accept the fact that the equivalence principle does not hold. If this is the case, a 
model with different mass parameters (violating the equivalence principle) would 
be more appropriate to describe the motion of test particles in the Solar system. 

The next step of our analysis of the phenomenological consequences of Matrix 
Gravity is to apply the kinematic model developed in [611 to the study of galactic 
rotations. It would be very interesting to understand if the flat rotation curves of 
galaxies can be explained without the concept of dark matter. We plan to investi- 
gate this question in a future work. 
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